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One of the most disastrous forms of collective human behaviour is 
the kind of crowd stampede induced by panic, often leading to fatal- 

ities as people are crushed or trampled. Sometimes this behaviour 

^" 

is triggered in life-threatening situations such as fires in crowded 

ON ■ 



buildings, ' at other times, stampedes can arise from the rush for 
seats 3,4 or seemingly without causes. Tragic examples within re- 
cent months include the panics in Harare, Zimbabwe, and at the 
Roskilde rock concert in Denmark. Although engineers are finding 
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ways to alleviate the scale of such disasters, their frequency seems 



to be increasing with the number and size of mass events. 2 ' 5 Yet, 
systematic studies of panic behaviour, 6 9 and quantitative theories 
capable of predicting such crowd dynamics, 51012 are rare. Here we 
show that simulations based on a model of pedestrian behaviour can 
provide valuable insights into the mechanisms of and preconditions 
for panic and jamming by incoordination. Our results suggest prac- 
tical ways of minimising the harmful consequences of such events 

1 
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and the existence of an optimal escape strategy, corresponding to a 
suitable mixture of individualistic and collective behaviour. 

Up to now, panics as a particular form of collective behaviour occuring in 
situations of scarce or dwindling resources 1,6 has been mainly studied from 
the perspective of social psychology. 7-9 Panicking individuals tend to show 
maladaptive and relentless mass behaviour like jamming and life-threatening 
overcrowding, 1-4 ' 8 which has often been attributed to social contagion 1,4 ' 8 
(see Brown 9 for an overview of theories). According to Mintz, 6 the observed 
jamming is a result of incoordination and depends on the reward structure. 

After having carefully studied the related socio-psychological literature, 6-9 
reports in the media and available video materials (see http : //angel . elte . 
hu/~panic/), empirical investigations, 1 ' 2,3 and engineering handbooks, 13,14 we 
can summarise the following characteristic features of escape panics: (i) Peo- 
ple move or try to move considerably faster than normal. 13 (ii) Individuals 
start pushing, and interactions among people become physical in nature, (hi) 
Moving and, in particular, passing of a bottleneck becomes incoordinated. 6 
(iv) At exits, arching and clogging are observed. 13 (v) Jams are building up. 7 
(vi) The physical interactions in the jammed crowd add up and cause danger- 
ous pressures up to 4,450 Newtons per meter, 2,5 which can bend steel barriers 
or tear down brick walls, (vii) Escape is further slowed down by fallen or in- 
jured people turning into "obstacles" . (viii) People show a tendency of mass 
behaviour, i.e., to do what other people do. 1 ' 8 (ix) Alternative exits are often 
overlooked or not efficiently used in escape situations. 1,2 

These observations have encouraged us to model the collective phenomenon of 
escape panic in the spirit of self-driven many-particle systems. Our computer 
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simulations of the crowd dynamics of pedestrians are based on a generalised 
force model, 15 which is particularly suited to describe the fatal build up of pres- 
sure observed during panics. 2-4 ' 5 We assume a mixture of socio-psychological 16 
and physical forces influencing the behaviour in a crowd: Each of N pedestri- 
ans % of mass rrii likes to move with a certain desired speed v® into a certain 
direction e°, and therefore tends to correspondingly adapt his or her actual 
velocity v t with a certain characteristic time Tj. Simultaneously, he or she 
tries to keep a velocity-dependent distance to other pedestrians j and walls 
W. This can be modelled by "interaction forces" fij and fiw, respectively 
In mathematical terms, the change of velocity in time t is then given by the 
acceleration equation 

dVi ^(t)e°(t) - Vj(t) ^ _ 

at T i j&i) w 

while the change of position Vi(t) is given by the velocity Vi(t) = dri/dt. 
We describe the psychological tendency of two pedestrians % and j to stay 
away from each other by a repulsive interaction force Ai exp[(r^- — dij)/Bj\ riy, 
where Ai and B{ are constants, dij = ||r*j — r*^ || denotes the distance between 
the pedestrians' centers of mass, and n^ = (n^-,n^-) = (t*j — Vj)/dij is the 
normalised vector pointing from pedestrian j to i. If their distance dij is 
smaller than the sum r^ = (r^ + Tj) of their radii r^ and Tj, the pedestrians 
touch each other. In this case, we assume two additional forces inspired by 
granular interactions, 17 ' 18 which are essential for understanding the particular 
effects in panicking crowds: a "body force" k{rij — dij) riij counteracting body 
compression and a "sliding friction force" K(rij — dij) Av^tij impeding relative 
tangential motion, if pedestrian i comes close to j. Herein, ty = (—n?j,n}j) 
means the tangential direction and At>*j = {vj — Vi) • tij the tangential velocity 
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difference, while k and k represent large constants. In summary, we have 
f i:j = {Ai exp[(r ij - dij)/Bi] + kg(r i:j - d i:j )} n*,- + Kg(r i:j - d^Av^ % , (2) 

where the function g(x) is zero, if the pedestrians do not touch each other 
(dij > rij), otherwise equal to the argument x. 

The interaction with the walls is treated analogously, i.e., if dm? means the 

distance to wall W, n iW denotes the direction perpendicular to it, and tiw 

the direction tangential to it, the corresponding interaction force with the wall 

reads 

f iW = {Ai exp[(rj - d iW )/Bi] + kg(n - d iW )} n iW - ng(ri - d iW ){vi ■ t iW ) t iW 

(3) 

Probably due to the fact that escape panics are unexpected and dangerous 
events, which also excludes real-life experiments, we could not find suitable 
data on escape panics to test our model quantitatively. This scarcity of data 
calls for reliable models. We have, therefore, specified the parameters as fol- 
lows: With a mass of rrii = 80 kg, we represent an average soccer fan. The 
desired velocity v® can reach more than 5 m/s (up to 10 m/s), 14 but the 
observed free velocities for leaving a room correspond to vf pa 0.6 m/s un- 
der relaxed, v® pa 1 m/s under normal, and v® ^1.5 m/s under nervous 
conditions. 13 A reasonable estimate for the acceleration time is Tj = 0.5 s. 
With Ai — 2 • 10 3 N and Bi = 0.08 m one can reflect the distance kept at 
normal desired velocities 14 and fit the measured flows through bottlenecks 14 , 
amounting to 0.73 persons per second for an effectively 1 meter wide door 
under conditions with v® pa 0.8 m/s. The parameters k = 1.2 • 10 5 kgs~ 2 and 
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k = 2.4 • 10 5 kgm _1 s _1 determine the obstruction effects in cases of physical 
interactions. Although, in reality, most parameters are varying individually, 
we chose identical values for all pedestrians to minimise the number of pa- 
rameters for reasons of calibration and robustness, and to exclude irregular 
outflows because of parameter variations. However, to avoid model artefacts 
(gridlocks by exactly balanced forces in symmetrical configurations), a small 
amount of irregularity of almost arbitrary kind is needed. This irregularity 
was introduced by uniformly distributed pedestrian diameters 2r^ in the inter- 
val [0.5 m, 0.7 m], approximating the distribution of shoulder widths of soccer 
fans. 

Based on the above model assumptions, we will now simulate several impor- 
tant phenomena of escape panic, which are insensitive to reasonable parameter 
variations, but fortunately become less pronounced for wider exits. 

1. Transition to Incoordination due to Clogging. The simulated outflow from 
a room is well-coordinated and regular, if the desired velocities v® = Vq are 
normal. However, for desired velocities above 1.5 m/s, i.e. for people in a rush, 
we find an irregular succession of arch-like blockings of the exit and avalanche- 
like bunches of leaving pedestrians, when the arches break (see Fig. la, b). 
This phenomenon is compatible with the empirical observations mentioned 
above and comparable to intermittent clogging found in granular flows through 
funnels or hoppers 17 ' 18 (although this has been attributed to static friction 
between particles without remote interactions, and the transition to clogging 
has been observed for small enough openings rather than for a variation of 
the driving force). 

2. u Faster-Is-Slower Effect" due to Impatience. Since clogging is connected 
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with delays, trying to move faster (i.e., increasing v®) can cause a smaller av- 
erage speed of leaving, if the friction parameter k is large enough (see Fig. lc, 
d). This "faster-is-slower effect" is particularly tragic in the presence of fires, 
where fleeing people reduce their own chances of survival. The related fatali- 
ties can be estimated by the number of pedestrians reached by the fire front 
(see http : //angel . elte . hu/~panic/) . 

Since our friction term has, on average, no deceleration effect in the crowd, 
if the walls are sufficiently remote, the arching underlying the clogging effect 
requires a combination of several effects: i) slowing down due to a bottle- 
neck such as a door and ii) strong inter-personal friction, which becomes 
dominant when pedestrians get too close to each other. Consequently, the 
danger of clogging can be minimised by avoiding bottlenecks in the construc- 
tion of stadia and public buildings. Notice, however, that jamming can also 
occur at widenings of escape routes! This surprising result is illustrated in 
Fig. 2. Improved outflows can be reached by columns placed asymmetrically 
in front of the exits, which also prevent the build up of fatal pressures (see 
http : //angel . elte . hu/~panic/). 

3. Mass Behaviour. We investigate a situation in which pedestrians are trying 
to leave a smoky room, but first have to find one of the invisible exits (see 
Fig. 3a). Each pedestrian % may either select an individual direction e* or 
follow the average direction (e°At))i of his neighbours j in a certain radius 
Ri, 19 or try a mixture of both. We assume that both options are weighted 
with some parameter pf. 

e°(t)=M\(l-p l )e l + Pt (e° 1 (t)} t } , (4) 
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where N(z) = z/\\z\\ denotes normalisation of a vector z. As a consequence, 
we have individualistic behaviour if p t is low, but herding behaviour if pi is 
high. Therefore, Pi reflects the degree of panics of individual %. 

Our model suggests that neither individualistic nor herding behaviour per- 
forms well (see Fig. 3b). Pure individualistic behaviour means that each 
pedestrian finds an exit only accidentally, while pure herding behaviour im- 
plies that the complete crowd is eventually moving into the same and probably 
blocked direction, so that available exits are not efficiently used, in agreement 
with observations. According to Figs. 3b and c, we expect optimal chances of 
survival for a certain mixture of individualistic and herding behaviour, where 
individualism allows some people to detect the exits and herding guarantees 
that successful solutions are imitated by the others. If pedestrians follow the 
walls instead of "reflecting" at them, we expect that herd following causes 
jamming and inefficient use of doors as well (see Fig. 1), while individualists 
moving in opposite directions obstruct each other. 

In summary, we have developed a continuous pedestrian model based on plau- 
sible interactions, which is, due to its simplicity, robust with respect to pa- 
rameter variations and suitable for drawing conclusions about the possible 
mechanisms beyond escape panic (regarding an increase of the desired veloc- 
ity, strong friction effects during physical interactions, and herding). After 
having calibrated the model parameters to available data on pedestrian flows, 
we managed to reproduce many observed phenomena including i) the build up 
of pressure, ii) clogging effects at bottlenecks, iii) jamming at widenings, iv) 
the "faster-is-slower effect", v) inefficient use of alternative exits (see Fig. 3d), 
and vi) initiation of panics by counterflows and impatience (i.e., seemingly 
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without any logical reason, see Fig. 1). Moreover, we are able to simulate 
situations of dwindling resources and estimate the casualties. Therefore, the 
model can be used to test buildings for their suitability in emergency sit- 
uations. It accounts for both, the different dynamics in normal and panic 
situations just by changing a single parameter p t = p. In addition, our sim- 
ulations suggest that the optimal behaviour in escape situations is a suitable 
mixture of individualistic and herding behaviour. 

We are now calling for complementary data and additional video material on 
escape panics to test our model quantitatively and compare it with alternative 
ones which, for example, include direction- and velocity-dependent interper- 
sonal interactions, specify the individual variation of parameters, study the 
effect of fluctuations, consider falling people, integrate acoustic information 
exchange, implement more complex strategies and interactions (also three- 
dimensional ones), or allow for switching of strategies. A superior theory 
would have to reproduce the empirical findings equally well with less parame- 
ters, reach a better quantitative agreement with data with the same number, 
or reflect additional observations. 
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FIG. 1. Simulation of pedestrians moving with identical desired velocity v® = vq 
towards the 1 m wide exit of a room of size 15 m x 15 m. a Snapshot of the scenario. 
Dynamic simulations are available at http://angel.elte.hu/~panic/. b Illus- 
tration of leaving times of pedestrians for various desired velocities vq. Irregular 
outflow due to clogging is observed for high desired velocities (vq > 1-5 m/s, see 
red plusses) . c Under conditions of normal walking, the time for 200 pedestrians to 
leave the room decreases with growing vq. Desired velocities higher than 1.5 m/s 
reduce the efficiency of leaving, which becomes particularly clear, when the outflow 
J is divided by the desired velocity (see d). This is due to pushing, which causes 
additional friction effects. Moreover, above a desired velocity of about vq = 5 m/s 
( — ), people are injured and become non-moving obstacles for others, if the sum of 
the magnitudes of the radial forces acting on them divided by their circumference 
exceeds a pressure of 1600 N/m. 5 

Due to the above "faster-is-slower effect", panics can be triggered by pedestrian 
counterflows, 2 which cause delays to the crowd intending to leave. This makes the 
stopped pedestrians impatient and pushy. One may describe this by increasing the 
desired velocity according to v®(t) = [1 — pi(t)]vf(0) +pi(t)vf iax , where v®(0) is 
the initial and vf 1 ^ the maximum desired velocity. The time-dependent parameter 
Pi(t) = 1 — Vi(t)/v^, where Vi(t) denotes the average speed into the desired direc- 
tion of motion, is a measure of impatience. Altogether, long waiting times increase 
the desired velocity, which can produce inefficient outflow. This further increases 
the waiting times, and so on, so that this tragic feedback can eventually trigger 
panics. It is, therefore, imperative, to have sufficiently wide exits and to prevent 
counterflows, when big crowds want to leave. 
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FIG. 2. Simulation of an escape route with a wider area (see also the Java 
applets supplied at http://angel.elte.hu/~panic/). a Illustration of the scenario 
with v® = Vo = 2 m/s. The corridor is 3 m wide and 15 m long, the length of the 
triangular pieces in the middle being 2x3 m = 6 m. Pedestrians enter the simulation 
area on the left-hand side with an inflow of J = 5.5 s _1 m _1 and flee towards the 
right-hand side, b Efficiency of leaving as a function of the angle </> characterising 
the width of the central zone, i.e., the difference from a linear corridor. The relative 
efficiency E = (Vi-e^}/vo measures the average velocity along the corridor compared 
to the desired velocity and lies between and 1 ( — ). While it is almost one (i.e., 
maximal) for a linear corridor (<fi = 0), the efficiency drops by about 20%, if the 
corridor contains a widening. This becomes comprehensible, if we take into account 
that the widening leads to disturbances by pedestrians, who expand in the wide area 
due to their repulsive interactions or try to overtake each other, and squeeze into the 
main stream again at the end of the widening. Hence, the right half of the illustrated 
corridor acts like a bottleneck and leads to jamming. The drop of efficiency E is 
even more pronounced, (i) in the area of the widening where pedestrian flow is most 
irregular ( — ), (ii) if the corridor is narrow, (iii) if the pedestrians have different or 
high desired velocities, and (iv) if the pedestrian density in the corridor is higher. 
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FIG. 3. Simulation of N = 90 pedestrians trying to escape a smoky room of area 
A = 15m x 15 m (grey) through two invisible doors of 1.5 m width, which have 
to be found with a mixture of individualistic and herding behaviour. Java applets 
are available at http://angel.elte.hu/~panic/. a Snapshot of the simulation 
with v ° = vq = 5 m/s. Initially, each pedestrian selects his or her desired walking 
direction randomly. Afterwards, a pedestrian's walking direction is influenced by the 
average direction of the neighbours within a radius of, for example, Ri = R = 5 m. 
The strength of this herding effect grows with increasing panic parameter Pi = p 
and increasing value of h = 7rR 2 p, where p = N/A denotes the pedestrian density. 
When reaching a boundary, the direction of a pedestrian is reflected. If one of the 
exits is closer than 2 m, the room is left, b Number of people who manage to escape 
within 30 s as a function of the panic parameter p. c Illustration of the time required 
by 80 individuals to leave the smoky room. If the exits are relatively narrow and the 
degree p of herding is small or large, leaving takes particularly long, so that only 
some of the people escape before being poisoned by smoke. Our results suggest 
that the best escape strategy is a certain compromise between following of others 
and an individualistic searching behaviour. This fits well into experimental data 
on the efficiency of group problem solving, 20 according to which groups normally 
perform better than individuals, but masses are inefficient in finding new solutions 
to complex problems, d Absolute difference \N\ — A^2 1 in the numbers Aq and N2 of 
persons leaving through the left exit or the right exit as a function of the herding 
parameter p. We find that pedestrians tend to jam up at one of the exits instead 
of equally using all available exits, if the panic parameter is large. 



